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Abstract 

The ground-state energy, the addition energies and the optical absorption spectra are derived 
for interacting polarons in parabolic quantum dots in three and two dimensions. A path integral 
formalism for identical particles is used in order to take into account the fermion statistics. The 
approach is applied to both closed-shell and open-shell systems of interacting polarons. Using a 
Jrj ■ generalization of the Jensen-Feynman variational principle, the ground-state energy of a confined 

iV-polaron system is analyzed as a function of N and of the electron-phonon coupling constant a. 
As distinct from the few-electron systems without the electron-phonon interaction, three types of 
spin polarization are possible for the ground state of the few-polaron systems: (i) a spin-polarized 
state, (ii) a state where the spin is determined by Hund's rule, (iii) a state with the minimal 
possible spin. A transition from a state fulfilling Hund's rule, to a spin-polarized state occurs when 
decreasing the electron density. In the strong-coupling limit, the system of interacting polarons 
turns into a state with the minimal possible spin. These transitions should be experimentally 
observable in the optical absorption spectra of quantum dots. 
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I. INTRODUCTION 



Many-electron states in quantum dots have been theoretically investigated by vari- 
ous approaches, e.g., the Hartree-Fock method 0, 0, El, the density-functional theory 
Ell the quantum Monte Carlo simulation the variational Monte Carlo 

method and the Pade approximation [l2[ a numerical diagonalization of the Hamil- 
tonian in a finite-dimensional basis The electron-phonon interaction was not taken 

into account in these investigations, although it can contribute significantly to both equilib- 
rium and non-equilibrium properties of quantum dots. For instance, the effects due to the 
electron-phonon interaction play a key role in the optical spectra of some quantum dots (see 
Ref . 2| and references therein) . Some characteristics revealed in the mid- infrared region of 
the experiment ally observed optical absorption spectra of high-T c cuprates [16ul7lll8l . ll9l . l20l ] 
were assigned |2l| to the polaron optical absorption at intermediate values of the electron- 
phonon coupling constant a. The infrared optical absorption band in the neodymium-cerium 
cuprate Nd 2 _ :r Ce a; Cu04_j / (NCCO), which was studied experimentally 2^] as a function of 
the electron density, has been associated with the polaron optical absorption. Also some 
peculiarities of the infrared optical absorption spectra of cuprates and ProNiO^w were in- 
terpreted in terms of a mixture of large and small polarons or bipolarons [23l . 

The theory of th e op tical conductivity of arbitrary-coupling single polarons has been 
developed in Refs. |25[ |26| (see also the review [ItJ and references therein) within the 
memory-function method based on the path- integral formalism |28|. Recently, the optical 
absorption of a gas of interacting polarons at weak coupling was investigated [i^, |3(| on 
the basis of a variational many-particle approach |3ll |. The optical absorption of interacting 
polarons in bulk semiconductors at arbitrary coupling strength has been treated in Ref. 
|32| within the random-phase a ppr oximation, using variational parameters obtained from 
Feynman's single-polaron model 28|| , but this treatment does not seem to be a self-consistent 
approach to the many-polaron problem. 

In contrast to the polaron mechanism of optical absorption in bulk (see Refs. (27I I33I ] 
and references therein), the polaron optical absorption in quantum dots, to the best of our 
knowledge, has not yet been widely studied. In order to investigate the ground state and 
the optical response of a system of interacting polarons in a quantum dot, it is crucial to 
take into account the fact that the system contains a finite number of identical particles 
(electrons or holes). Indeed, the thermodynamic properties of systems with a finite number 
of identical particles might substantially deviate from those obtained within the grand- 
canonical formalism (see, e.g., [3|3^). The variational path-integral method for identical 
particles 0, 3^, 3^, 37J provides a useful tool for investigating interacting quantum many- 
body systems with a fixed (few or many) number of N particles. An outline of the method 
is given in Ref. [HJ , where we sketched the calculation of the ground-state energy of a fixed 
number of interacting polarons, which form a closed-shell system in a quantum dot. 

In the present paper the calculation of the ground state and of the optical conductivity is 
performed for both closed-shell and open-shell systems of interacting polarons in a quantum 
dot. In Section El we derive an upper bound to the free energy of a finite number of 
interacting polarons confined in a parabolic quantum dot in three dimensions and in two 
dimensions. In Section HTT1 we discuss the numerical results for the ground-state energy and 
for the addition energy of this system. In Section HVl the optical conductivity of interacting 
polarons in a quantum dot is derived on the basis of the memory-function method. The 
numerical results for the optical conductivity are discussed in Section El The last section 
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IVII contains conclusions. 



II. THE PARTITION FUNCTION AND THE FREE ENERGY OF A MANY- 
POLARON SYSTEM 

A. Interacting polarons in a quantum dot 

We consider a system of N electrons with mutual Coulomb repulsion and interacting 
with the lattice vibrations. The system is assumed to be confined by a parabolic potential 
characterized by the frequency parameter Qq. The total number of electrons is represented 
as N = N a , where N a is the number of electrons with the spin projection a = ±1/2. The 
electron 3D (2D) coordinates are denoted by Xj )CT with j = 1, • • • , N a . The bulk phonons 
(characterized by 3D wave vectors q and freq uencies u; q ) are described by the complex 
coordinates Q q , which possess the property |39| 

Q q = Q-q- (1) 

The full set of the electron and phonon coordinates are denoted by x = {xj j(T } and Q = {Q q } • 
Throughout the sections UT1 and HTT1 the Euclidean time variable r = it is used, where t is 
the real time variable. In this representation the Lagrangian of the system is 

L (x, Q; x, 0) = L e (x, x) - V c (x) + L ph (Q, q) + L e _ ph (x, Q) , (2) 
where L e (x, x) is the Lagrangian of an electron with band mass mj in a quantum dot: 



<7=±l/2 j = l 

Vc (x) is the potential energy of the electron-electron Coulomb repulsion in the medium 
with the high-frequency dielectric constant : 

M*) = E EEf, — z — r ^ 

,/=±l/2i=l 1=1 Z£ °° |X ^ CT Xi < CT/| 

-^pft ^0; 0*; Qi Q*^j is the Lagrangian of free phonons: 

L ph (Q,Q) =-\Y&Q* + ^QkQ^ Q^§- ( 5 ) 

q 

Further, L e _ p h (x, Q,Q*) is the Lagrangian of the electron-phonon interaction: 

L e - ph (X, Q)=-J2 {j^) 1 W-qPq^ (6) 
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where p q is the Fourier transform of the electron density operator: 

p*= Y, » ,,x "'- ( 7 ) 

a=±l/2 j=l 

V q is the amplitude of the electron-phonon interaction. In this paper, we only consider 
electrons interacting with the long-wavelength longitudinal optical (LO) phonons with a 
dispersionless frequency cj q = co> LO , for which the amplitude V q is |26| 

_ ^ LO / 2V2na \ ^ ( h \ 1/4 

where a is the electron-phonon coupling constant and V is the volume of the crystal. 

We treat a canonical ensemble, where the numbers N a are fixed. The partition function 
Z ({N a } , 0) of the system can be expressed as a path integral over the electron and phonon 
coordinates: 

Z({JV„},0) = E ^gg^ / dx f CS(r) JdQj%Q(T) e -*VWt , (9) 

where S 1 [x (r) , Q (r)] is the "action" functional: 

1 / - -\ 

S[x{t),Q{t)]=-j:J L(k,Q;x,Q)dT. (10) 

The parameter (3 = 1/ (kgT) is inversely proportional to the temperature T. In order to take 
the Fermi-Dirac statistics into account, the integral over the electron paths {x(r)} in Eq. 
(JHJ) contains a sum over all permutations P of the electrons with the same spin projection, 
and £p denotes the parity of a permutation P. 

The action functional (P"U|) is quadratic in the phonon coordinates Q. Therefore , th e path 
integral over the phonon variables in Z ({N a } , 0) can be calculated analytically [39|. As a 
result, the partition function of the electron-phonon system © factorizes into a product 

z(w,^vw,fln iM J tWi) <") 

of the partition function of free phonons with a partition function Z p ({N a } , 0) of interacting 
polarons, which is a path integral over the electron coordinates only: 



^ IRW)= ?/wfe/ fi f w<l)t " 4W ' (12) 



The functional 

-ft/3 



i r np 

S P [x (r)] = -- J [L e (x (r) , x (r)) + V c (x (r))] rfr 



ft/3 ft/3 

Y^J dT J dT sinh(^ LO /2) * (T) P -« (T } (13) 

q 
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results from the elimination of the phonon coordinates and contains the "influence phase" of 
the phonons (the last term in the right-hand side). It describes the phonon-induced retarded 
interaction between the electrons, including the retarded self-interaction of each electron. 
The free energy of a system of interacting polarons F p ({N a } , (3) is related to their partition 
function (J12p by the equation: 

F p ({N a } ,P) = -~hiZ p ({N a } , (3) . (14) 

At present no method is known to calculate the non-gaussian path integral (I12|) analyti- 
cally. For distinguishable particles, the Jensen-Feynman variational principle |39| provides a 
convenient approximation technique. It yields a lower bound to the partition function, and 
hence an upper bound to the free energy. 

The formulation of a variational principle for the free energy for a system of identical par- 
ticles is a non-trivial problem. However, it can be shown j^f that the path-integral approach 
to the many-body problem for a fixed number of identical particles can be formulated as a 
Feynman-Kac functional on a state space for N indistinguishable particles, by imposing an 
ordering on the configuration space and by the introduction of a set of boundary conditions 
at the boundaries of this state space. The path integral (in the imaginary-time variable) 
for identical particles was shown to be positive within this state space. This implies that a 
many-body extension of the Jensen-Feynman inequality was found, which can be used for 
interacting identical particles (Ref.j36[, p. 4476). A more detailed analysis of this varia- 
tional principle for both local and retarded interactions can be found in Ref. [3?}. It is 
required that the potentials are symmetric with respect to all permutations of the particle 
positions, and that both the exact propagator and the model propagator are antisymmetric 
(for fermions) with respect to permutations of any two electrons at any moment in time. 
This means that these propagators have to be defined on the same configuration space. 
Keeping in mind these constraints, the variational inequality for identical particles takes the 
same form as the Jensen-Feynman variational principle: 



F P < F + I (S p - S ) So , (15) 



where So is a model action with corresponding free energy F . The angular brackets mean 
a weighted average over the paths 



((•))s - rtfp . r p x _ „ r , y • ( 16 ) 



B. Model system 

We consider a model system consisting of N electrons with coordinates x = {x J)(J } and 
Nf fictitious particles with coordinates y = {y^} in a harmonic confinement potential with 
elastic interparticle interactions as studied in Ref. [38]. The Lagrangian of this model system 
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takes the form 



N a 2 N ° N °' 



L M (i, y; x, y) = ^ ^ (x 2 ff + fi 2 x 2 ff ) + ^ - 

a j = l <r,cr' j = l 1=1 

- f E W + "M) - 1 E E E - y) 2 • < 17 > 

j = l (7 j = l i=l 

The frequencies Q, to, Qf, the mass of a fictitious particle mj, and the force constant k 
are variational parameters. Clearly, this Lagrangian is symmetric with respect to electron 
ermutations. Performing the path integral over the coordinates of the fictitious particles 
the partition function Z ({N a } , (3) of the model system of interacting polarons becomes 
a path integral over the electron coordinates: 



Z„({iV4,/?) = £^£-! / fi^flxMe-^M, (18) 



r - '1/2^-1/2 

with the action functional S$ [x (r)] given by 

-hp N, 



s » I* Ml = h / EE? W. M + « 2 4, M] *■ 

1 /"^ ^ ^ m^ 2 2 



cr.cr' i=l (=1 

hp hp 


where X is the center-of-mass coordinate of the electrons, 

x =^£5>>- ( 2 °) 

cr j = l 

The details of the analytical calculation of the model partition function ([18)1 are described 
in Appendix A. 

After substituting the model action functional (|19|) into the right-hand side of the varia- 
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tional inequality (fT5j). we obtain an upper bound to the free energy F v , 
F var ({N a },(3) 



N 

F ({N a } ,/3) + ^(nl-n 2 + n. 2 ) ( J>J (o) 



So 



^ <* (")>* + £ ^ fa (1, 01 W, /*)-"] 

k 2 N 2 N f f f cosh[tt f (\T-T'\-hp/2)} 



+ ^mr,} dT J dT ' s mh( W ) ^ m^w* 



l^ql /* j /" , / Cosh[a; L o (\t-t'\- hp/2)] /|f , n . , 01 , 

7 rf V dr sinh(^ LO /2) 9(d,T-T\{N a },(]). (21) 



o 



Here, g (q, r — t'\ {N a } , (3) is the two-point correlation function for the electron density 
operators: 

^(q,r|{iV (T },/?) = (p q (r)p_ q (0)) 5o . (22) 

Both the free energy and the correlation functions of the model system can be calculated 
analytically using the generating function technique j23|. In the zero-temperature limit 
(P — > oo) , the variational free energy (|2"Tj) becomes an upper bound E° ar ({N a }) to the 
ground-state energy E° of the system of interacting polarons. The details of the calculation 
of the correlation functions are given in the Appendix B. 



III. GROUND-STATE ENERGY AND ADDITION ENERGY OF INTERACTING 
POLARONS 

For the numerical calculations, we use effective atomic units, where h, the electron band 
mass rrib and e/y/e^ have the numerical value of 1. This means that the unit of length 
is the effective Bohr radius a* B = f^E^j (m^e 2 ), while the unit of energy is the effective 
Hartree[I3| H* = mi > e A / (ffe 2 ^). These units allow to present results for quantum dots with 
and without the electron-phonon interaction on the same scale. Therefore, for confined 
polarons they are more convenient than the usual polaron units, where the unit of length is 
a p = [h/ (m^Lo)] 1 ^ 2 , and the energy is measured in units of the LO-phonon energy /j^lo 
In terms of the dimensionless parameters a and rj = Eoc/eq, where is the static dielectric 
constant, the following relations exist between both systems of units: 

ap = V2« H* fa P \ 2 _ 2a 2 

a* B 1-v hu LO \a* B J ~ (I-?]) 2 ' 1 J 

In Fig. 1, the total spin 5* of a system of interaction polarons in their ground state is 
plotted as a function of the number of electrons in a 3D quantum dot for different values of the 
confinement frequency Qq, of the electron-phonon coupling constant a and of the parameter 
rj. As distinct from few-electron systems without the electron-phonon interaction, three types 
of spin polarization are possible for the ground state, which should be distinguishable from 
each other using, e.g., capacity measurements. 
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(i) Except for the strong-coupling case and for the low-density case, the filling in the 
ground state is as follows: in an open shell, with less-than-half filling, each new electron is 
added with one and the same spin, so that the total spin (in the shell under consideration) 
is maximal in accordance with Hund's rule jZoj . As soon as half- filling is achieved with 
electrons possessing a certain spin, each new electron is added with the spin opposite to 
that in the group of electrons providing the aforementioned half-filling. When the number 
of electrons corresponds to the number of states in the shell under consideration, the shell 
becomes closed, and the total spin is zero. This mode of filling is referred to as Hund's rule 
for a quantum dot. Hund's rule means, that the electrons in a partly filled upper shell build 
up a minimal possible number of pairs in order to minimize the electron-electron repulsion. 
For a quantum dot with hflo = 0.5 H* at a = and at a = 0.5, the shell filling always 
obeys Hund's rule, as shown in Fig. 1(a). 

(ii) With decreasing confinement frequency Qq at a given number of electrons, the electron 
density lowers. For densities smaller than a certain value, it can happen that the ground 
state is a state with a maximal total (in all shells) spin. In this state, the electrons are filling 
consecutively all the single-electron states with one and the same spin and are referred to 
as spin-polarized. The examples are the states at a = for N = 4 and N = 10 in Fig. 
1(6). A spin-polarize d gr ound state precedes the formation of a Wigner crystal when further 



lowering the density |14l . 141 



(iii) In the strong-coupling case ( a > 1 and r] <C 1), it can happen that the ground 
state is a state with a minimal total spin (0 for even number of electrons and ~ for odd 
number of electrons). This is the case when -due to the phonon- mediated electron-electron 
attraction- pairing of electrons with opposite spins occurs, analogous to a singlet bipolaron 
ground state in bulk. The examples are the states at a = 5 and rj — 0.1 for N in the range 
from 4 to 6 in Fig. 1(a) and for N in the range from 4 to 10 in Fig. 1(b). This trend to 
minimize the total spin is a consequence of the electron-phonon interaction, presumably due 
to the fact that the phonon-mediated electron-electron attraction overcomes the Coulomb 
repulsion. With an increasing number of electrons, at a certain value of N, such states with 
a minimal total spin cease to form the ground state, and the shell filling abruptly returns 
to that prescribed by Hund's rule [see a jump in the spin at a = 5 when N changes from 10 
to 11 in Fig. 1(b)]. This jump is analogous to a transition from states with paired electrons 
(like superconducting states) to another type of states with unpaired electrons (like normal 
states). 

The addition energy A (N), which is the variation of the chemical potential when putting 
an extra electron into a quantum dot, is defined as 0,0 

A (iV) = E° (N + 1) - 2E° (N) + E° (N - 1) . (24) 

Fig. 2 presents the addition energy in a 3D quantum dot as a function of the number 
of electrons. The structure of A (N) clearly manifests the shell structure of a quantum 
dot. The most pronounced peaks in the addition energy occur for closed-shell systems with 
N = 2,8, 20. The peaks in A (N) at iV = 5 and N = 14 obtained within the present 
approach for relatively weak electron-phonon coupling correspond to the systems with the 
half-filled upper shell (see Fig. 2(a) for a = and a = 0.5). In these cases the total 
spin for the upper shell takes its maximal possible value, in accordance with Hund's rule. 
At sufficiently large values of a, the electron-phonon interaction substantially modifies the 
addition energy. In the strong-coupling case, the peaks corresponding to half-filled shells 
become less pronounced, while those corresponding to closed-shell systems become more 
prominent as compared to the weak-coupling case. 
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To the best of our knowledge, the addition energy for parabolic quantum dots was ob- 
tained using the density functional theory (DFT) (see, e.g., Refs.0, IE H|) only without the 
electron-phonon interaction. Our results for the addition energy for a 3D quantum dot as 
a function of N in the particular case a = are very close for iV $C 12 to those calculated 
within the DFT Q| with an optimized effective potential and a self-interaction correction 

Hi. 



The panels a and b in Fig. 3 represent, respectively, the total spin and the addition 
energy for interacting polarons in a 2D parabolic GaAs quantum dot with the confinement 
parameter hflo = 0.5 H* « 7.67 meV. The pronounced peaks in A (N) at N = 2, 6, 12, 20, . . . 
correspond to the closed-shell systems, for which the total spin equals zero. In accordance 
with Hund's rule, the upper shell is filled in such a way that the total spin of electrons in 
this shell takes the maximal possible value. Therefore for the half-filled upper shell (at N = 
4, 9, 16, . . .) maxima of the total spin occur as a function of N. At these electron numbers, 
the addition energy manifests peaks, which are less pronounced than those corresponding 
to closed shells. 

The inset to Fig. 3(b ) shows the experimental data for the addition energy in a cylindric 
GaAs quantum dot [43]. As seen from Fig. 3, the peak positions for the addition energy of 
interacting polarons in a 2D parabolic quantum dot agree well with the experimental results 
for the addition energies of cylindrical quantum dots. The height of the calculated peaks of 
the addition energy falls down as the shell number increases, which is qualitatively consistent 
with the experimentally observed behavior |43j]. The peaks in A (N) corresponding to the 
half-filled shells are weaker than those for the closed shells both in the experiment 43[ and 
in our theory. 

IV. OPTICAL CONDUCTIVITY 

For a system of interacting polarons in a parabolic confinement potential, we calculate 
the real part of the optical conductivity within the memory-function app roach. For a single 
polaron at arbitrary coupling strength it was developed in Ref s. | 25l l2q . For a polaron gas 
in the weak-coupling limit, this technique was applied in Ref. |44j |. 

In the present paper we extend the memory-function approach to a system of arbitrary- 
coupling interacting polarons. Since the optical conductivity relates the current J (t) per 
electron to a time-dependent uniform electric field E (t) in the framework of linear response 
theory, we have to return to the real-time representation in the path integrals. The Fourier 
components of the electric field are denoted by E w : 



1 



OC 



E(t) = — / E u e-**du, (25) 

and the similar denotations are used for other time- dependent quantities. The electric 
current per electron J (t) is related to the mean electron coordinate response R (t) by 

JW = -^. (26) 

and hence 

J w = ieioK^. (27) 
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Within the linear-response theory, both the electric current and the coordinate response are 
proportional to E w : 

^>E„„ (28) 



where a (u) is the conductivity per electron. Because we treat an isotropic electron-phonon 
system, a (u) is a scalar function. It is determined from the time evolution of the center-of- 
mass coordinate: 



*(*)^((f>i(') 



(29) 

ki=l ll s 

The symbol {{{*))) s denotes an average in the real-time representation for a system with 
action functional S: 



«(.)»«, = dx dxo dx' / Dx(t) / Dx'(t)e 



iS[x(t),x'(t)] / \ / 



(x |p (t )| x c 



to— oo 



x„ 



(30) 

where (x \p (to) I x o) is the density matrix before the onset of the electric field in the infinite 
past (to — > — oo). The corresponding action functional is pH Hil] 



S [x (t) , x' (t)] = / [L e (x (t) , x (t) , t) - L e fx (t) , x' (t) , t) 1 dt' - ih$ [x (t) , x' (t)] , (31) 



where L e (x, x, t) is the Lagrangian of N interacting electrons in a time-dependent uniform 
electric field E (t) 



JV CT 



L e (x,x,t) = EE 

(7 jf = l 



The influence phase of the phonons (see, e.g., Ref. |46j|) 



No- N a , 



- ex ,>- E (*) -EEE 



ct,o-' i=i z=i 



(32) 



2 S 

$[x( S ),x'( S )] = -^^L da ' [p q (s) - // q («)] 



-oo — oo 



x \T:(s-s>) Pci (s>)-T^(s-s>)p'(s') 



(33) 

describes both a retarded interaction between different electrons and a retarded self- 
interaction of each electron due to the elimination of the phonon coordinates. This functional 
contains the free-phonon Green's function: 



Auit 



-iuit 



+ 



(34) 



The equation of motion for R (t) can be derived by analogy with that described in Ref. j!?} : 

d 2 R (t) 



m b - 



dt 2 



+ m b n 2 Q R (t) + eE (t) = F ph (t) 



(35) 
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where F p h (t) is the average force due to the electron-phonon interaction, 



F ph (t) = -ReJ2 [ ds T: lq (t - s) «p q (t) p_ q (s)» s . (36) 

q 

The two-point correlation function ((p q (t) p_ q (s)))s should be calculated from Eq. (JBTIJ) 
using the exact action (J3T|) . but like for the free energy above, this path integral cannot 
be calculated analytically. Instead, we perform an approximate calculation, replacing the 
two-point correlation function in Eq. (j3Hj) by ((p q (t) p_ q (s)))s Q , where So [x (t) , x' (t)] is 
the action functional with the optimal values of the variational parameters for the model 
system considered in the previous section in the presence of the electric field E(t). The 
functional S [x (t) , x' (t)} is quadratic and describes a system of coupled harmonic oscillators 
in the uniform electric field E (t). This field enters the term — eE (t) • ^ CT YlfZi x j,°- m the 
Lagrangian, which only affects the center-of-mass coordinate. Hence, a shift of variables to 
the frame of reference with the origin at the center of mass 



(37) 



x n (t) =x n (t) + R(t) 

X ^ (t)=a n (t) + R(f) 

results in [46] 

((p q (t)P- q m So = ((pAt)P-As))) s J E=o e i(l - m) - Kis)] - (38) 

This result (I38j) is valid for any quadratic model action So. 

The applicability of the parabolic approximation for N = 1 is confirmed by the fact, 
that for the polaron ground-state energy, the results of the Feynman approach |28| are very 
close to the values obtained using other reliable methods j3, EH 5^, 57, Hlf. Thus, a self- 



induced polaronic potential, created by the polarization cloud around an electron, is rather 
well described by a parabolic potential whose parameters are determined by a variational 
method. For N = 2, the lowest known values of the bipolaron ground-state energy are 
provided by the path-integral variational method with parabolic potentials both in bulk 
|59| and for confined systems |60|, |6l( for realistic values of a. The aforesaid approximation 
for the right-hand side of Eq. (J36|) is a direct generalization of the all-coupling approach 
(2H . liH r53l | to a many-polaron system. For weak coupling, our variational method is at 
least of the same accuracy as the perturbation theory, which results from our approach at a 
special choice of the variational parameters. For strong coupling, an interplay of the electron- 
phonon interaction and the Coulomb correlations within a confinement potential can lead 



to the assemblage of polarons in multi-polaron systems. As shown in Refs. [60|, |61| for a 
system with N = 2, the presence of a confinement potential strongly favors the bipolaron 
formation. Our choice of the model variational system is reasonable because of this trend, 
apparently occurring in a many-polaron system with arbitrary iV for a finite confinement 
strength. 

The correlation function (p q (t) p_ q (s)) 5o | E corresponds to the model system in the ab- 
sence of an electric field. For t > s, this function is related to the imaginary-time correlation 
function g (q, t\ {N a } , 0) , described in the previous section: 

«Pq (*) P- q (s))) So | E=0 t>s = 9 (q, i (t ~ s) I {iV CT } , 0) . (39) 
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Using the transformation (|37p and the relation (|39|). one readily obtains 



F ph (t) = -ReJ2 / T: LO (t - s) ^""Wlj (q, • (t - a) | {iV CT } , /3) <fe. (40) 

Within the framework of the linear-response theory, the external electric field E (t) is a 
small perturbation, so that R (t) is a linear functional of [E (t')]L'<t • Expanding the function 
e iq-[R(t)-R(s)] j n right-hand side of Eq. (|4Uj) in powers of [R(£) — R(s)] up to the first- 
order term, we obtain the Fourier component F p h (u) of the force due to the electron-phonon 
interaction which is proportional to R w . As a result, the optical conductivity can be expressed 
in terms of the memory function x ( w ) (cf. Refs. [25, 26]), 



Re a (u) = o, 41 

m 6 [ w a_ng_Hex(w)] 3 + [Imx(w)] a 

where x ( u ) is given by 

2 °° 

= £ / dt " x) Im [ T: - (t) 9 (q ' ul {N ° } ' ^ • (42) 

<J 6 

It is worth noting that the optical conductivity (}4*T]) differs from that for a translationally 
invariant polaron system both by the explicit form of x ( w ) an d by the presence of the term 
in the denominator. For a — > 0, the optical conductivity tends to a 5-like peak at uj = Qq, 

7TP 2 

lim Re a (u) = 5 (u - fi ) • (43) 

a->o 2?m 



For a translationally invariant system fi ~^ 0, and this weak-coupling expression 
reproduces the "central peak" of the polaron optical conductivity 

In the zero-temperature limit, the memory function of Eq. (|42j) is derived in the analytical 
form of Eq. (jB.15|) in Appendix B for 3D and 2D interacting polarons. 



V. RESULTS ON THE OPTICAL CONDUCTIVITY 



Due to the confinement, the electron motion in a quantum dot is fully quantized. Hence, 
when a photon is absorbed, the electron recoil can be transferred only by discrete quanta. As 
a result, the optical conductivity spectrum of a system of interacting polarons in a quantum 
dot is a series of <5-like peaks as distinct from the optical conductivity spectrum of a bulk 
polaron (25I . |26| . These peaks are related to the internal polaron excitations. 

Because Im x (w) = for all frequencies except for a discrete set of combinatorial fre- 
quencies (IB.16|) . the peaks in the optical conductivity (}4*Tj) are positioned at the frequencies 
which are given by the roots of the equation 

u 2 -n 2 -Re X (u;) = 0, (44) 
which are denoted as ( f2i, f2 2 > • • • ) • 
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One of these roots is close to the variational parameter , which is the eigenfrequency of 
the motion of a polaron as a whole. It satisfies the inequality Q2 < ^0 because the polaron 
effective mass is larger than that of a bare electron. Q2 is close to Qq in the weak-coupling 
case and decreases with increasing a. Hence, it tends to zero in the limit Qq — *> 0. The peak 
in Re a (u) corresponding to this root can be considered as the zero-phonon line, which is 
an analogue of the "central peak" of the polaron optical conductivity (jjl . 26l |. The peaks 



of Recr (c<j) determined by the other roots of Eq. can be attributed to transitions into 
excited states of the many-polaron system. 

The changes of the shell filling schemes, which occur when varying the confinement fre- 
quency, manifest themselves in the spectra of the optical conductivity. In Fig. 4, optical 
conductivity spectra for N = 20 polarons are presented for a quantum dot with the param- 
eters of CdSe: a = 0.46, rj = 0.656 [49j and with different values of the confinement energy 
MV In this case, the spin-polarized ground state changes to the ground state satisfying 
Hund's rule with increasing htt in the interval 0.0421if* < fiCl < 0.0A22H*. 

In the inset to Fig. 4, the first frequency moment of the optical conductivity 

_ J Q °°ujRea(uj)duj 
J Re cr [00) duo 

as a function of MIq shows a discontinuity, at the value of the confinement energy corre- 
sponding to the change of the shell filling schemes from the spin-polarized ground state to 
the ground state obeying Hund's rule. This discontinuity should be observable in optical 
measurements. 

In Fig. 5, the first frequency moment (|45|) is plotted as a function of the number of 
electrons for a CdSe quantum dot with fl = 0.143co>lo (corresponding to hfl ~ 0.04 H*). 
The total spin of the system as a function of N is shown in the inset. As a general trend, (u) 
decreases with increasing N, with kinks corresponding to the ground-state transitions from 
states obeying Hund's rule with iV = 3, 9, and 18, into spin-polarized states with N = 4, 
10, and 19, respectively. 

In Fig. 6, optical conductivity spectra are plotted for several values of the confinement 
frequency for N = 10 polarons in a quantum dot with a = 2, rj = 0.6. These values of a 
and 7] are typical for the high-T c superconducting cuprates of the NCCO family 29]. In 
the "weak- confinement" region (f2 = 0.6u;lo an d — 0.8^0) the zero-phonon peak is 
expressively dominant over the other peaks. 

When the confinement frequency parameter passes through the value fl =_u; L o, the so- 
called "confinement- phonon resonance" [3^ occurs. In this case, the peaks at Q^, k = 1,2,3 
have comparable oscillator strengths. The position Q2 of the second peak is substantially 
shifted from the LO phonon frequency culo Moreover, the intensities of the phonon-assisted 
transitions increase as compared to the "weak- confinement" case. This resonance has a clear 
analogy with the magneto-phonon resonance (see, e.g., Ref. {50}), as far as the energy levels 
of an electron in a parabolic confinement are similar to the Landau levels of an electron in 
a magnetic field. 

With further increasing f2 , when Q > uj^q, the dominant part of the optical conductivity 
spectrum shifts to higher frequencies. For instance, at Qq = 1.4ig>lo the most intensive peak 
is that with Cl 3 . The intensities of the peaks, beginning with the second peak, increase 
in comparison with the intensities of their "weak-confinement" analogs. The positions of 
the zero-phonon line and the subsequent peaks are substantially shifted from the "weak- 
confinement" values towards higher frequencies. These effects are a manifestation of the 
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mixing of the zero-phonon state with different excited states of the many-polaron system. 
A similar behavior of the optical absorption spectra at and above the magneto-phonon 
resonance is explained by the mixing of zero-phonon and one-phonon quantum states [Hof . 

The shell structure for a system of interacting polarons in a quantum dot is clearly 
revealed when analyzing the addition energy and the first frequency moment of the optical 
conductivity in parallel. In Figs. 7 and 8, we show both the function 

e(N)= (u)\ N+1 -2(u)\ N +{co)\ N _ 1 , (46) 

and the addition energy A (N) for interacting polarons in different 3D quantum dots. 

As seen from Figs. 7 and 8 for quantum dots of CdSe and with a = 3j5l|, respectively, 
distinct peaks appear in (N) and A (N) at the "magic numbers" corresponding to closed- 
shell configurations at iV = 8, 20 for the states obeying Hund's rule in panels a, b and to 
half-filled- shell configurations at iV = 10, 20 for the spin-polarized states in panels c, d of 
Fig. 8. In the case when the shell filling scheme is the same for different N (see panels a, b 
in Figs. 7, 8, where the filling obeys Hund's rule), each of the peaks of (N) corresponds 
to a peak of the addition energy. In the case when the shell filling scheme changes with 
varying iV (panels c, d in Figs. 7, 8), the function (N) exhibits pronounced minima for N 
corresponding to the change of the filling scheme from the states, obeying Hund's rule, to 
the spin-polarized states. 

It follows that measurements of the addition energy and the first frequency moment of the 
optical absorption as a function of the number of polarons in a quantum dot can reflect the 
difference between open-shell and closed-shell configurations. In particular, the closed-shell 
configurations may be revealed through peaks in the function (JV). The filling patterns 
for a many-polaron system in a quantum dot can be determined from the analysis of the 
first moment of the optical absorption for different numbers of polarons. The appearance of 
minima in the function (N) will then indicate a transition from the states which are filled 
according to Hund's rule to the spin-polarized states. 



VI. CONCLUSIONS 

We presented a formalism for calculating the ground-state energy and the optical conduc- 
tivity spectra of a system of N interacting polarons in a parabolic confinement potential for 
arbitrary electron-phonon coupling strength. The path integral treatment of the quantum 
statistics of indistinguishable particles 0, allows us to find an upper bound [3g to the 
ground-state energy of a finite number of polarons. The parameters from the variational 
procedure are used as input for the calculation of the optical conductivity spectrum of the 
system. 

Two types of transitions were found for N polarons confined in a parabolic potential, 
with the corresponding ground states characterized by different values of the total spin. 
In the weak-confinement regime, the polaron system is in the spin-polarized state. When 
increasing the confinement frequency Qq, the system goes into a state obeying Hund's rule 
at a specific value of Q . For a strongly coupled system of interacting polarons, a third type 
of state appears, for which the total spin takes its minimal value. The analysis is performed 
for both closed-shell and open-shell systems. 

The calculations of the optical conductivity spectra for a finite number of polarons in a 
quantum dot are based on the memory-function approach. The dependence of the optical 
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conductivity spectra on the confinement parameter f2o reveals a resonant behavior for Qo ~ 
c^lo Transitions between states with different values of the total spin manifest themselves 
through discontinuous changes of the optical conductivity spectra and of the addition energy 
function of the number of electrons. 
The first frequency moment of the optical conductivity as a function of the number of 
electrons clearly shows the transition between the spin-polarized ground state of interacting 
polarons in a quantum dot and the ground state obeying Hund's rule, and it also can be used 
to discriminate between open-shell and closed-shell configurations. Optical measurements 
are therefore suggested as possible tools for examining the shell structure of a system of 
interacting polarons. 
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A. PARTITION FUNCTION OF THE MODEL SYSTEM 

In this appendix we discuss the analytical calculation of the partition function 
Zq ({^o-} > 0) [Eq. (HBJ)] f° r the model system of interacting polarons. It can be expressed in 
terms of the partition function Zm , Nf, (3) of the model system of interacting electrons 

and fictitious particles with the Lagrangian L M [Eq. (fT7j) ] as follows: 

MW,.^™, (A.1) 
z f(Nf,w f ,P) 

where Zf (Nf,Wf,f3) is the partition function of fictitious particles, 



Z f (N f , (3) = j^, (A.2) 



with 



w f = JQj + kN/rrif (A.3) 



and D=3(2) for 3D(2D) systems. The partition function Zm {{N a } , Nf, (3) is the path 
integral for both the electrons and the fictitious particles: 



/f PyL f f y 

dx D5l(t) / dy / Dy (r) e- SM[ ^ T) ^ T)] (A.4) 
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with the "action" functional 

-h/3 



1 f p 

S M [x (r) , y (r)] = -- J L M (x, y; x, y) dr, (A.5) 

where the Lagrangian is given by Eq. (|17|) . 

Let us consider an auxiliary "ghost" subsystem with the Lagrangian 

L 9 (± g , Y g , X g , Y 9 ) = -?f (XJ + w^ g ) - (yJ + w)Y]) (A.6) 
with two frequencies w and Wf, where w is given by 

w = sjtt 2 - Nu 2 + kN f /m b . (A.7) 
The partition function Z g of this subsystem 

Z g = j dX g I dY 9 J DX g (r) J DY 9 (r) exp {-S g [X g (r) , Y g (r)]} , (A.8) 

with the "action" functional 

h/3 

S g [K g (r) , Y g (r)] = ~\ f L a (x g , X fl , Y„ Y,) rfr (A.9) 

o 

is readily calculated: 

Z * = 7 777777777 777— 77J- (A. 10) 



[2sinh(^)] 



2s inh(^ 



The product Z g Z M of the two partition functions Z g and Z M {{N a } , Nf, (3) is a path 
integral in the state space of N electrons, Nf fictitious particles and two "ghost" particles 
with the coordinate vectors X g and Y g . The Lagrangian Lm of this system is a sum of Lm 
and L n 



L M (x, y, X g ,Y g ; x, y, X g , Y fl ) = L M (x, y; x, y) + L 9 (x ff , Y ff , X ff , Y ff ) . (A. 11) 

The "ghost" subsystem is introduced because the center-of-mass coordinates in L M can 
be explicitly separated more easily than in Lm- This separation is realized by the linear 
transformation of coordinates, 

x j> = x j> + X - X s , , . 

v = v' +Y-Y 

where X and Y are the center-of-mass coordinate vectors of the electrons and of the fictitious 
particles, correspondingly: 

X =4££^ Y =wX y - (A - 13) 

a j = l ■> j = l 
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Before the transformation ([A. 12)1 . the independent variables are (x, y, X , Y s ) , with the 
center-of-mass coordinates X and Y determined by Eq. (jA.13|) . After this transformation 
the independent variables can be considered to be (x',y',X, Y), where the coordinates 
(Xg,Y g ) obey the equations 



a 3=1 ' 3=1 



(A. 14) 



A substitution of Eq. (|A.13J) into Eq. (|A.11|) results in the following 3 terms: 

L M (x', y', X, Y; x', y', X, y) = L w (x', x') + L Wf (y', y') + L c (X, X; Y, y) , (A. 15) 

where L w fx , x'J and L Wf fy', y'J are Lagrangians of non-interacting identical oscillators 
with the frequencies w and Wf, respectively, 



er=±l/2 3=1 



AT, 



(A.16) 
(A.17) 



The Lagrangian Lc (^X, X; Y, YJ describes the combined motion of the centers-of-mass of 
the electrons and of the fictitious particles, 



L C (X,X;Y,Y 



mbN 'x 2 + fi 2 X 2 



m f N f ( Y 2 + w 2 f Y 2 )+kNN f X-Y, (A.18) 



with 



n = \ Kl 2 + kN f /m b . 



(A.19) 



It is reduced to a diagonal quadratic form in the coordinates and the velocities by a standard 
transformation for two interacting oscillators: 



X 
Y 



1 



y/m b N 
1 



(air + a 2 R) , 



-a 2 r + a x R) 



(A.20) 



with the coefficients 



ai 



X 



"i + x" 


1/2 


[i-xl 




, a-2 = 


2 


2 



1/2 



fi 2 -n 2 



ft 2 -ft 2 ) +4 7 2 



1/2' T — ^1 



AA 



/ 



rribmf 



(A.21) 
(A.22) 
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The eigenfrequencies of the center-of-mass subsystem are then given by the expression 




(A.23) 



As a result, four independent frequencies f^i, Q 2 , w and Wf appear in the problem. Three of 
them (Qi, Q 2 , w) are the eigenfrequencies of the model system. Q\ is the frequency of the 
relative motion of the center of mass of the electrons with respect to the center of mass of 
the fictitious particles; Q 2 is the frequency related to the center of mass of the model system 
as a whole; w is the frequency of the relative motion of the electrons with respect to their 
center of mass. The parameter w / is an analog of the second variational parameter w of the 
one-polaron Feynman model. Further, the Lagrangian (jA.18|) takes the form 



(A.24) 



leading to the partition function corresponding to the combined motion of the centers-of- 
mass of the electrons and of the fictitious particles 



1 



j c 



[ 2s inh(»)] D [2sinh(»)] 



D ■ 



(A.25) 



Taking into account Eqs. (jA.lOjl and (|A.25|) . we obtain finally the partition function of 
the model system for interacting polarons 



Z o ({N a },0 



sinh (^) sinh 



-, D 



sinh (») sinh (£fa) 



Z F ({N a },w,P). 



Here 



Z F ({N a } ,w,0) = Z F (N 1/2 , w, 13) Z F (iV_ 1/2 , w, 0) 



(A.26) 



(A.27) 



is the partition function of = N1/2+N-.U2 non-interacting fermions in a parabolic confine- 
ment potential with the frequency w. The analytical expressions for the partition function 
of N a spin-polarized fermions (N a , w, 0) were derived in Ref. 34 . 



B. TWO-POINT CORRELATION FUNCTION. MEMORY FUNCTION. 



The two-point correlation function (|2~2~|) is represented as the following path integral: 



E 



^o(W^)V%^-i/2! 



JdxJ Dx(r)e- 5o ^ T )]p q (r)p- q (0). (B.l) 
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We observe that g (q, r\ , (3) can be rewritten as an average within the model "action" 
Sm [x (t) , y (r)] of interacting electrons and fictitious particles: 



g(q,T\{N a },fl) 



Z M ({N a },N f ,(3)^N 1/2 \N_ 1/2 \ 



E 



Px 



x y dx y D5c(r) I dy / Dy(r)e 

X Pq (T) P-q (0) • 



-5 M [x(r),y(r)] 



(B.2) 



Indeed, one readily derives that the elimination of the fictitious particles in (jB.2j) leads 
to (jB.lj) . The representation (jB.2j) allows one to calculate the correlation function 
g (q, r | {Ncr} , P) in a much simpler way than through Eq. (jB.ljl . using the separation of 
the coordinates of the centers of mass of the electrons and of the fictitious particles. This 
separation is performed for the two-point correlation function (jB.2j) by the same method as 
it has been done for the partition function (|A.4j) . As a result, one obtains 



g(q,T\{N a },P)=g( q ,r\{N CT },(3) 



(exp[*q- (X(r)-X(a))]) 



Sc 



(exp [iq- (X g (r) - X 9 (a 



(B.3) 



where g (q, t\ {N a } , 0) is the time-dependent correlation function of N non-interacting elec- 
trons in a parabolic confinement potential with the frequency w, 

^(q,r|{iV (T },/3) = (p q (r)p_ q (0)) 5ro . (B.4) 

The action functional S w [x T ] is related to the Lagrangian L w (x, x) [Eq. (IAT61) ] 

hp 

S w [x r ] = ^Jl w (x,x) dr. (B.5) 



The averages in (jB.3jl are calculated using Feynman's method of generating functions 

sinh 



(exp [iq- (X (r) - X (<r))]) Sc = exp 



hq 2 



Nm 



2 sinh 

2 



2 



2 



(exp [iq- (X g (r) - X g (<t))]) 5o = exp 



^ S inh(«) 
hq 2 sinh (^) sinh 



Nrrih 



w 



sinh (^) 



The two-point correlation functi on q (q, r| {Ao-} , /3) is derived using the generating-function 
technique for identical particles [34|. After the path integration, the following expression is 



obtained: 



g(q, -ir\{N a },(3) 



n,o\n' ,cr' 



+ E K eiq - X )nJ 



exp 



x [A (n, <r| {iV ff } ,/3)-f 2 (n, a; n', cr| {N a } 



(B.6) 
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where (e iqx ) nn , is the one-electron matrix element, 

(e lQ X ) nn , = / e iq ' x C (x) (x) dx. (B.7) 

For a 3D quantum dot, ip a (x) is the eigenfunction of a 3D oscillator with the frequency w 
[see, e.g., Ref. The index n denotes the set n = (n, /, m) , where n is the number of the 

energy level e n = hw (n + 3/2), / is the quantum number of the orbital angular momentum 
and m is the quantum number of the z projection of the orbital angular momentum. Simi- 
larly, for a 2D quantum dot, ip n (x) is the eigenfunction of a 2D oscillator with the frequency 
w. 

The one-electron distribution function f\ (n,a\N a , /3) is the average number of electrons 
with the spin projection a at the n-th energy level, while the two-electron distribution 
function f 2 (n, a; n', a'\ {N a } , (3) is the average product of the numbers of electrons with the 
spin projections a and a' at the levels n and n' . These functions are expressed through the 
following integrals (cf. Ref. [35^): 



fx (n, a\N a , (3) = - - * / / (e n , 9) $ (9, f3, N a ) d9, (B.8) 



2ttZ f (N ff ,w,P) 

— TT 

J f (e n , 0) f (e n r, 9) $ (9, (3, N a ) d9, if a' = a; 

B.9) 



h (n, a; ri, a'\ {N a } ,/3)={ 

h (n, a\N a , [3) h (n, a'\N a ,,f3) , if a' ± a 



with the notations 

^(9,(3,N a )=exp 



In (1 + e*-*-*") ~N a {£ + i9) 



n=0 



(B.IO) 



f ^ s -s^rwri - (B11) 

The function / (e, 6 1 ) formally coincides with the Fermi-Dirac distribution function of the 
energy e with the "chemical potential" (£ + z#) / /3. 

From here on we consider the zero-temperature limit, for which the integrals (jB.8|) and 
()B.9|) can be calculated analytically. The result for the one-electron distribution function is 




h (n^lAiV.)!^ = m 0, n > ( B .12) 



9£ 



According to ()B.12|) . L a is the number of the lowest open shell, and 

\(n + 1) {n + 2) (3D) , 



n + 1 (2D) . 



is the degeneracy of the n-th shell. iV^ is the number of electrons in all the closed shells 
with the spin projection a, 



La 1 



N - _ \ \L a (L a + 1) (L a + 2) (3D), 

^ = ^ " 1 6 I L(T (L„ + 1) (2D) . ( R13 ) 



n=0 
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The two-electron distribution function / 2 (n, a; n', a'\ {N a } ,f3) at T = takes the form 



/ 2 (n,a;n>'|/?,{^})l^oo 

A (n, a% NJl^ h (ri, a'\(3 } N^. 

1, 
0, 

N-N LtT N N l 1 



9L a 



►oo > n ^ n' or (j ^ d' 

a = a' and n = n' < L a \ 
a = a' and n = n' > L a \ 
a = a' and n = n' = L„. 



(B.14) 



Finally, using the two-point correlation function Eq. (jB.7|) . the one-electron ()B.12|) and 
the two-electron (|B.14|) distribution functions, the memory function of Eq. (|4*2*|) can be 
represented in the unified form for 3D and 2D interacting polarons: 



e^+o 3nN V 4 J \ A J 



x 



oo oo oo 

Pl=0p2=0 P3=0 
oo oo 



-1) 



pMW \nq 1 a 



pi 



NQ 2 A 



I>2 



1 \ 



NwA 



X 



X 



E E E ^ ^ ff i w ' $ - ^ ( n > a ' m ' a i w> ' 



m=0 n=0 cr 



w-ulo — [pifh +P2^2+(P3 -m+ra)to]+ie w+^lo+Pi^i +P2^2+(P3 -m+n)w+ie 



-V 



v Wlo+Pi^i+P2^2+(P3— m+n)w 
\n—m+l+k 



) 



*E E 

(=0 k=n— m+l 

n + D-1 

n — k 



k r{ Pl +p 2 + P3 + k + i + l) 



km 



i 

wA 



l+k 



X 



+ 



2k 

k — Z — n + m 
l 



i 



w-w L 0-(pi^l+P2n2+_P3'"')+ie w+o; L o+pini+p2^2+P3«'+i£ 



oo oo 



X 



X 



l^LO+Pl^l+P2^2+P3™ , 

i) fe+ 'r( Pl +p 2 +p 3 + fc + / + |) 



m=0 n=0 CT,fr' 



EE 

fe=0 1=0 

n + D-1 
n — k 



kin 

m + D - 1 

m — Z 



1 



(B.15) 



where D = 2, 3 is the dimensionality of the space, V denotes the principal value, A is defined 
as A = [X)i=i + (N — 1) /w] /AT, fi 1; f2 2 , and w are the eigenfrequencies of the model 

system, a± and a 2 are the coefficients of the canonical transformation which diagonalizes the 
model Lagrangian (JT7j) derived in Appendix A. 

The typical spectra of the real and imaginary parts of the memory function \ (iS) , are 
plotted in Figs. 9(a) and 9(b), respectively. According to Eq. (jB.15|) . the poles of Rex (k>) 
and the 5-like peaks of [- Imx(w)] are positioned at the combinatorial frequencies % m , 
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which are linear combinations of the LO-phonon frequency and three eigenfrequencies 

ujkim = ^lo + k0 1 + 10 2 + mu>, (B.16) 

with integer coefficients k, I, m = 0, 1, Each combinatorial frequency uJkim corresponds 

to a phonon-assisted transition to an excited state of the model system. 

The roots of the equation (jUJ), which provide the peaks in the optical conductivity, 

\0±, O2, ■ ■ .V are indicated in Fig. 9 a by the vertical arrows. For the chosen parameters, the 

peak at 0\ is the zero-phonon line. Fig. 9a also reveals peaks of Kea(uj) with frequencies 
in between two neighboring discrete values of uj n hi : e.g., at 2 . Following the physical 
interpretation of the memory function in Refs. [25|, |2y, |53[, these peaks can be related to 
transitions into excited states of the many-polaron system. 
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Figure captions 

Fig. 1. Total spin of a system of interacting polarons in a parabolic quantum dot as a 
function of the number of electrons for HQq = 0.5H*(a) and for HQq = 0.1H*(b). 

Fig. 2. Addition energy of a system of interacting polarons in a parabolic quantum dot 
as a function of the number of electrons for hQo — 0.5H*(a) and for HQq — 0.1H*(b). 

Fig. 3. The total spin (a) and the addition energy (b) of a system of interacting polarons 
in a 2D parabolic GaAs quantum dot as a function of the number of electrons for HQq = 
0.5H*. Inset: the experimentally observed addition energy vs number of electrons in a 
cylindrical GaAs quantum dot for two values of the diameter D = 0.5/xm and D = 0.44/im 

Fig. 4. Optical conductivity spectra of JV = 20 interacting polarons in CdSe quantum 
dots with a = 0.46, r] = 0.656 for different confinement energies close to the transition 
from a spin-polarized ground state to a ground state obeying Hund's rule. Inset: the first 
frequency moment (u) of the optical conductivity as a function of the confinement energy. 

Fig. 5. The first frequency moment (u>)of the optical conductivity as a function of the 
number of electrons for systems of interacting polarons in CdSe quantum dots with a = 0.46, 
r] = 0.656 and 0.143o;lo (M^o ~ 0.04H*). Open squares denote the spin-polarized ground 
state; full dots denote the ground state, obeying Hund's rule; open triangles denote the 
ground state of the third type, with more than one partly filled shells, which is not totally 
spin-polarized. Inset: the total spin of the system of interacting polarons as a function of 
N. 

Fig. 6. Optical conductivity spectra of JV = 10 interacting polarons in a quantum dot 
with a = 2, 1] = 0.6 for several values of the confinement frequency from Qq = 0.6u;lo to 
Qo = 1.4u;lo- The spectrum for Qq = c^lo corresponds to the confinement-phonon resonance. 

Fig. 7. The function O (JV) and the addition energy A (N) for systems of interacting 
polarons in CdSe quantum dots with a = 0.46, r\ = 0.656 for hfl = 0.1H* (panels a, b) and 
for Q = 0.04H* (panels c,d). Open squares denote the spin-polarized ground state; full 
dots denote the ground state, obeying Hund's rule; open triangles denote the ground state 
of the third type, with more than one partly filled shells, which is not totally spin-polarized. 

Fig. 8. The function O (JV) and the addition energy A (JV) for systems of interacting 
polarons in quantum dots with a = 3, rj = 0.25 and Qo = uj^q (panels a, b) and with a = 3, 
?7 = 0.3 and Q = 0.5o;lo (panels c, d). 

Fig. 9. Real (a) and imaginary (b) parts of the memory function x (^) f° r a system 
of interacting polarons in a quantum dot for JV = 10, a = 2, 77 = 0.6, and fl = 0.6^0- 
The dashed line in panel a represents (uj 2 — Qq) ju). The vertical arrows in panel a indicate 
the roots of Eq. ([44)1 . The height of peaks in panel b represents the relative intensity of the 
5-like peaks of [— Im% (u) /uj}. 
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